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We present the results of a lattice QCD calculation of the pseudoscalar meson decay constants 
fK, fo and fo,, performed with Nf = 2 dynamical fermions. The simulation is carried out with 
the tree-level improved Symanzik gauge action and with the twisted mass fermionic action at 
maximal twist. With respect to our previous study [|lj], here we have analysed data at three val- 
ues of the lattice spacing (a ~ 0.10fm.0.09fm,0.07fm) and performed the continuum limit, and 
we have included at a = 0.09 fm data with a lighter quark mass (mj^ ~ 260MeV) and a larger 
volume (L ~ 2.7fm), thus having at each lattice spacing L > 2.4fm and m^^L > 3.6. Our result 
for the kaon decay constant is /k = (157.5 ± 0.8|i(oi. ±3.31 iwr.)MeV and for the ratio Jk/ fn = 
1.205 ±0.006|,,„r.± 0.0251 syst., in good agreement with the other Nf ~ 2 and Nf — 2+1 lattice 
calculations. For the Z) and meson decay constants we obtain//) = (205 ±7|s,a,. ±7|.„sr.)MeV, 
in good agreement with the CLEO-c experimental measurement and with other recent Nf = 2 and 
Nf = 2+1 lattice calculations, and fu, = (248 ±l>\stat. ± 8|i„r.)MeV that, instead, is 2.3a be- 
low the CLEO-c/BABAR experimental average, confirming the present tension between lattice 
calculations and experimental measurements. 
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1. Introduction 

We present a lattice QCD determination [Q] of the pseudoscalar meson decay constants fx, 
fo and fo,, performed at three values of the lattice spacing a ~ O.lOfm, 0.09 fm, 0.07 fm (corre- 
sponding to j8 = 3.8,3.9,4.05), with pion masses down to nijt ~ 260MeV and volumes such that 
L > 2.4 fm and nijiL > 3.6, at each lattice spacing.' The continuum extrapolation turns out to be 
crucial for an accurate determination of fo and fo, since cutoff effects induced by the charm mass 
(~ &{a^ml) ~ 0.04^0.09) are not small. 

The calculation is based on the gauge field configurations generated by the European Twisted 
Mass Collaboration (ETMC) with the tree-level improved Symanzik gauge action and the twisted 
mass action at maximal twist, discussed in detail in ref. [Q], with the Nf = 2 dynamical quarks taken 
to be degenerate in mass. The use of the twisted mass fermions in the present calculation turns out 
to be beneficial, since the pseudoscalar meson masses and decay constants, which represent the 
basic ingredients of the calculation, are automatically improved at ^{a) and the determination 
of pseudoscalar decay constants does not require the introduction of any renormalization constant. 

In order to investigate the properties of the K, D and mesons, we simulate the sea and 
valence light (u / d) quark mass in the range O.lmf^'- <Hi< O.Smf^'-, where nis^^'^' is the physical 
strange mass, the valence strange quark mass within 0.9 nif''-^*' < /^^ < 1. 5 and the valence 

charm quark mass within O.Smf'"^' < /^c < l.Snic''^''''' , being ruc'^-'^'" the physical charm mass. 

The statistical accuracy of the meson correlators is improved by using the so-called "one-end" 
stochastic method, implemented in ref. [Q], which includes all spatial sources. Statistical errors on 
the meson masses and decay constants are evaluated using the jacknife procedure and statistical 
errors on the fit results which are based on data obtained from different configuration ensembles 
are evaluated using a bootstrap procedure. 

In the present analysis we study the dependence of the pseudoscalar decay constants on the 
meson masses instead of quark masses, as in our previous analysis [p at fixed value of the lattice 
spacing. The study in terms of meson masses is simpler here, where data at different values of 
the lattice spacing are involved, since it does not require the introduction of the quark mass renor- 
malisation constant (Z„, = Zp that would be necessary if the analysis were carried out in terms 
of physical quark masses. The dependence of the decay constants on the meson masses is studied 
simultaneously with the dependence on the lattice spacing, through a combined fit where terms 
of ff{a^) and ff{a^ jiq) are added to the functional forms predicted by Chiral Perturbation Theory 
(ChPT). In studying the kaon and the Dg meson sectors, we have treated the strange quark mass 
either as a light quark by using SU(3)-ChPT for kaons and SU(3)-Heavy Meson ChPT (HMChPT) 
for Ds mesons, or by treating only the light u/d quarks as light. In the latter case we use SU(2)- 
ChPT and the interpolation to the physical strange quark is performed linearly. This is justified, 
since our simulated values of the strange quark mass are quite close to the physical strange mass. 



2. Determination of the kaon decay constant 
Fit based on SU(3)-ChPT 

'The results of a lattice calculation of the vector and tensor decay constants in the kaon sector, performed with 
Nf = 2 twisted mass QCD by a subgroup of our ETM Collaboration, have been recently presented in ref. 



2 



fK, fo ond fo^ from Nf = 2 tniQCD 



Cecilia Tarantino 



0,43 
0,4? - 

0,41 - 

0,40 - 
t 0,39 - 
0,38 - 
0,37 - 
0,36 - 
0.35 - 



P^3.8, L/a^24 
P-3.9, L/a-24 
P=3.9, L/a=32 
P=4.05, L/a=32 
SU(3)-ChFr 
SU(2)-ChPT 



0.2 0,3 0,4 0.5 0.6 0.7 0.8 



1.0 



Figure 1: Lattice results for r^fK = '"o//'s(M;,M/jMi) as a function of the simulated pion meson mass square 
(rom^i)^ = {''oMpsiHi , Hi , Hi))^ ■ We display data with /i^ fixed to the simulated mass that corresponds to a 
ss meson of mass roMps{lJ.i,IJ.s,IJ.s) = 1.63. The dashed (dotted) curves represent the SU(3)- (SU(2)-) ChPT 
extrapolation to the physical pion mass, both at fixed lattice spacing (upper curves) and in the continuum 
limit (lower curve). The shift of the physical results, square (SU(3)) and diamond (SU(2)) dots, below the 
continuum limit curves comes from the interpolation to the physical strange quark mass. 

We perform a combined fit of tiie data available for the pseudoscalar decay constants at the three 
values of the lattice spacing, by using the functional form predicted by continuum NLO SU(3)- 
ChPT with the addition of an tff{a^) and ^{a^pLs) term. The expansion for fps{pisea,l^y!a'i:l^^^i), 
where jXsea and jX^^^ denote generically the sea and valence quark masses respectively, reads 



rofpsil^ul^hl^s) 



+biit,u+bss^ss) , 



^ ^11 In ^11 - ^ In ^ss - In 2(§/, 



(2.1) 



with the variables i^'s expressed as a function of meson masses as = M^^{pii,^i,^j)/ {AnfY . 
The parameter / is one of the low energy constants (LECs) entering the chiral Lagrangian at the LO, 
bii and bss are related to the NLO LECs It, and /4, whereas the coefficients A and Am parameterize 
discretization effects. The data for the pion decay constants fps{li.i, 1X1,1X1), also included in the 
analysis, are fitted through the same functional form with jXs = jXi. The fit is performed in 
units of the Sommer parameter ro [Q], for which values of ro/a at the three lattice spacings have 
been extracted in ref. [Q] from the analysis of the static potential. The values of ro and of the u/d 
quark mass in the continuum limit are extracted here by using as experimental input the pion mass 
and decay constant. 

We have also tried a fit with NNLO terms proportional to ^^,^ss ^ii^ss^ to take into account 
higher order chiral corrections, but their coefficients are found to be compatible with zero, showing 
that the NLO formula is satisfactory. It is worth noting that, in general, NLO SU(3)-ChPT doesn't 
describe well the data for the pseudoscalar decay constant up to the kaon sector |||, Indeed, in our 
previous analysis [|l|] performed in terms of quark masses we found that NNLO terms were needed 
to fit the data. The difference here is that the analysis is performed in terms of meson masses, 
where the replacement t, = Mp^/ {^Tif)^ effectively resums higher order chiral contributions and 
we find that the fit based on NLO SU(3)-ChPT is accurate enough to describe the pseudoscalar 



decay constant up to the kaon sector [10]. 
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As far as discretization effects are concerned, they are found to be at the level of ~ 2% on the 
finest lattice. The quality of the combined chiral and continuum extrapolation of the kaon decay 
constant is illustrated in fig. |l|. 
Fit based on SU(2)-ChPT 

As recently pointed out in ||8|], a good accuracy of NLO SU(3)-ChPT is in general not guaranteed 
in the kaon sector, where a safer approach consists in avoiding the chiral expansion in terms of the 
strange quark mass and applying therefore SU(2)-ChPT. In particular, the SU(2)-ChPT formula for 
the decay constant can be obtained from the SU(3)-ChPT expression by expanding in the quark 
mass ratio /i/ / pLs and including the strange quark mass dependence in the SU(2) LECs. In order to 
have small values for /i/ / pLs that justify the use of SU(2)-ChPT, we consider in the analysis only 
data with pLs in the range of the physical strange quark mass and we exclude, at each value of j8, 
data with values of /X/ corresponding to iriji > 500MeV . 

In order to perform a combined chiral and continuum extrapolation fit, we introduce the terms 
taking into account the discretization effects in the fitting formula, that reads 



l-^^„ln5,+Z.W5, 



(2.2) 



where the SU(2) LECs /(^' and b^^'^ are functions of the strange quark mass. We can safely assume 
a linear dependence on the strange mass, since the simulated /i, masses are close to the physical 
strange quark mass. 

From fig. | it is evident that SU(3)- and SU(2)-ChPT fits yield very similar results, f^^^^ = 
(157.9 ± 0.8) MeV and /^^^^^ = (157.1 ±0.7)MeV, with similar statistical uncertainty, and with a 
bit smaller chi-squared from SU(2)-ChPT {x^/d.o.f. = 0.8) than from SU(3)-ChPT {x'^/d.o.f. = 
1.2). 

Results for Jk and fK/fn 

The analyses based on SU(3)- and SU(2)-ChPT provide results that are in very good agreement and 
that have similar statistical uncertainties. We choose to average them and to quote their deviation 
from the average as the systematic uncertainty due to the chiral extrapolation. In order to estimate 
the uncertainty coming from discretization effects, we consider the difference between the values 
taken by the kaon decay constant after performing the continuum limit and at the finest lattice 
spacing, i.e. a ~ 0.07fm for j8 = 4.05. This latter value reads /^|4.o5 = 160.8MeV and turns 
out to be ~ 2% above the continuum limit result. In the present analysis, finite size effects (FSE) 
are estimated by using NLO ChPT [ 1 1 1 and are found to be negligible in the kaon sector. They 
only affect the determination of ro from the pion decay constant. As expected in simulations with 
L > 2.4fm and mj^L > 3.6, the FSE in the pion sector are found to be well under control, as 
confirmed by the comparison of data available at two volumes (L = 2.0fm and L = 2.7 fm, with 
a = 0.09 fm) and by the compatibility between the results for ro determined here by treating FSE 



within NLO ChPT Ql 1[ ] and those obtained in [|4j] by using the resummed formulae of ref. [|L 
Concerning the uncertainty due to the quenching of the dynamical strange quark, we believe that 
such an effect is smaller than the other systematic uncertainties estimated above, as suggested by 
the good agreement between recent Nf = 2 and Nf = 2 + 1 lattice determinations [^. 
We thus quote our final results for the kaon decay constant and the ratio /k/Zk 

h = (157.5 ±0.8 1, s,,,. ±0.4|,;„-, ±3.3Ui,,,)MeV = (157.5 ±0.8|,,«,. ±3.3|,,,t.)MeV, 
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fK/h = 1.205 ±0.006|,s,«,. ±0.003|,ft,>. ±0.025| 



discr. 



1.205 ±0.006L,„, ±0.0251 



syst. 5 



(2.3) 



where the total systematic uncertainty is obtained by adding in quadrature the chiral and discretiza- 
tion errors. Our result for the ratio fK/ fn turns out to be in good agreement with most of other 
recent lattice determinations [^] based on simulations with Nf = 2 and N f = 2 + \ dynamical 
fermions and with the average based on the determination of V^s from /T^s decays [|T^. 

3. Determination of the D and decay constants 



In order to determine the D and Dg meson decay constants we essentially proceed as for the 
kaon sector, that is we analyse simultaneously data at various values of the lattice spacing, by 
performing for the pseudoscalar decay constants combined fits of the meson mass dependence and 
discretization terms, in units of r^. The simulated values jlc of the charm quark mass are close to the 
physical charm quark mass (O.Swic''-^*' < jlc < 1.5m''''-^'* ), so that the interpolation to the physical 
value represents a very safe step. Moreover, at those large values of the meson masses the FSE are 
completely negligible. On the other hand, the discretization effects depending on the charm mass 
have to be taken into account in the fit, being parameterically of ff{a^ jj.^) ~ 0.04^-0.09. 

The functional forms describing the mass dependence of the decay constants, used to fit the 
data in the D and sectors are those predicted by the HMChPT where, as in the case of the kaon 
sector, one can consider two different approaches to treat the strange quark based on SU(2)- and 
SU(3)-ChPT, respectively. 

We determine fo and fo, by introducing two ratios that exhibit a smoother chiral behaviour: 



(3.1) 



Ik Ik fD^/mE 

In addition, discretization effects in the ratio R2 vanish in the limit of exact SU(3) symmetry. In 
the ratios the decay constants are multiplied by the square roots of the meson masses to reconstruct 
the quantities that are finite in the infinite mass limit. The Heavy Quark Effective Theory (HQET), 
in fact, predicts for a Heavy(//)-light(/) meson: fniVWii = A + B/Mhi + G(\ /Mjj^). Though the 
charm quark mass is far from the infinite mass Umit, in our analysis we can safely assume such a 
dependence for the D mesons, where only a small interpolation to the physical charm quark mass 
is needed. Moreover, since the contribution of the 1/Mhi correction is small, we consider the 
dependence on the light meson masses only in the leading term, by using the HMChPT prediction. 
Fit based on SU(3)-HMChPT 

Within the analysis based on SU(3)-HMChPT, the functional forms used to fit the ratios Ri and R2 



of eq. ( |3.lD read 



rn Ri 



:Ci- ( l+Cs^+Ce-M,, 



'0 



a2 



'3 1 
+ In ^1, + ^i. In ^is + ^ ^1! In 



l + 3g 



2^,M^,s + -{^n 



■2^,.)lni§.« 



+ 



Ri 



■Ci- 



2 2 
1+Cii-+Ci2-(M, 



'0 



-M,i) 



l + (l + 3# 



roMc, 
.2. /3 



(3.2) 



(^//ln<^//-(^/,.ln<^/,. --(<^// 



-2,^,,)ln,^,,) + 



■ ^ ^,1 In ^,1 + 5, In ^i, + i ^i, In ) + Cg ^// + C9 i,. 



+ 



Cio 



ro \M, 



1 



Mci, 
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Figure 2: Lattice results for the ratio Ri defined in the text, as a function of ttie pion mass square m\ = 
MpsiUhUhUi)^, in units ofro. We display data with jXs and jXe fixed to the simulated masses that correspond 
to a ss and a Ds meson of masses rQMps{jj.i,jj.s,IJ's) — 1-63 and r^Mcs = roMps{lJ.i,IJ.s,IJ.c) — 4.41. The dashed 
(dotted) curves represent the SU(3)- (SU(2)-) HMChPT extrapolation to the physical pion mass, both at fixed 
lattice spacing (upper curves) and in the continuum limit (lower curve). The shift of the physical results 
below the continuum limit curves comes from the interpolation to the physical strange and charm masses. 

where the coefficients C1-C12 are free fit parameters. The HMChPT parameter g cannot be deter- 
mined from the fit, which is almost insensitive to it, and is thus constrained to g = 0.5 [14]. By 
varying g in the range [0.3,0.7] we find that the fit results vary within the statistical uncertainty. 

The dependence on the pion mass square for the two ratios Ri and R2 is softer than for the 
quantities fo.^JmD, and /dV^^d shown for Ri in fig. ^. 

Fit based on SU(2)-HMChPT 

We obtain the SU(2)-HMChPT formulae for the decay constants and consequently for the ratios 
Ri and R2 from the SU(3)-HMChPT expressions by expanding in jj^i / jj.^ and including the strange 
quark mass dependence in the SU(2) LECs. They read 



D 



D 



1 • 



2 2 

\+Ai^+A2—Mc, 
4 ro 



(3.3) 



R2 = D4[l+A,^+A4-iMc,-M,,) 1 



ia2 



-^iiln^ii+Ds^ii + 



1 



1 



where the coefficients Di, D2, D4 and D5 are fitted as linear functions of the strange quark mass. 

From fig. |it is evident that the SU(3)- and SU(2)-HMChPT fits yield similar results (/^^'^^ = 
(205±7)MeV, /^^'^^ = (204±7)MeV and /^^^^^ = (249±3)MeV, /^^'^^ = (247±3)MeV) 
with practically the same statistical uncertainty and quality of the fit. 
Results for fp, fp^ and fp JfD 

The analyses based on SU(3)- and SU(2)-HMChPT provide results for fp, fp^ and for the ratio 
/dJ/d that are in perfect agreement, with very similar statistical uncertainties. As for the kaon 
decay constant, we choose to average the SU(3)- and SU(2)-HMChPT results and to quote their 
deviation from the average as the systematic uncertainty due to the chiral extrapolation. In order 
to quote the uncertainty coming from discretization effects, we consider the difference between 
the values assumed by fp, fp^ and fpjfp after the continuum limit has been taken and at the 
finest lattice spacing. The latter values read /z)|4.05 = 212MeV, /dJ4.05 = 256MeV, that are ~ 3% 
above the continuum limit results and (/d,//d)|4.05 = 1-207 that is 0.3% below the value in the 
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continuum limit. FSE effects are invisible in the D and sectors and can be neglected in the 
estimate of the systematic uncertainties. Concerning the uncertainty due to the quenching of the 
dynamical strange quark, as for the kaon sector, we believe that such an effect is smaller than the 
other systematic uncertainties conservatively estimated above. 

We therefore quote as our final results for the D and decay constants and the ratio foj fu 

fo = (205 ± 7 1 ,tat. ± 1 1 chir. ± lUiscr. ) McV = (205 ± 7 1 ,n,t. ± 7 1 ) MeV , 

/z), = (248±3|,,„,. ± ±8U„,,)MeV = (248 ± 3|,,„,. ± 8|,,,0MeV , (3.4) 

foJfD = 1.211 ± 0.0351,,,,. ± 0.003 1,;,,-, ±0.004 U,,,, = 1.211 ±0.035|,,,,. ±0.005|,,,,. . 



The result obtained for fo is in very good agreement with the CLEO-c measurement 015|], fj^ 



rexp. 



(205.8 ± 8.5 ±2.5) MeV and with other A^j = 2 and A^j = 2 ± 1 lattice calculations Even more 
interesting is the comparison of our f^^ result with the experiments and to other lattice results. 
The average of the recent CLEO-c JItI ] and BABAR [18] experimental measurements is f'j^'^' = 



(277 ± 9) MeV, tipically higher than the values indicated by lattice calculations and with a possible 
explanation as an effect of New Physics []T9|]. Our result is in the bulk of other recent lattice 
determinations [ p^ ] and confirms at the 2.3a level the tension with the experimental average. 
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